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Microscopic approach of fission dynamics 
applied to fragment kinetic energy 
and mass distributions in 238 U 
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The collective dynamics of low energy fission in 238 U is described within a time-dependent for- 
malism based on the Gaussian Overlap Approximation of the time-dependent Generator Coordi- 
nate Method. The intrinsic deformed configurations of the nucleus are determined from the self- 
consistent Hartree-Fock-Bogoliubov procedure employing the effective force D1S with constraints on 
the quadrupole and octupole moments. Fragment kinetic energy and mass distributions are calcu- 
lated and compared with experimental evaluations. The effect of the collective dynamics along the 
fission paths and the influence of initial conditions on these distributions are analyzed and discussed. 

PACS numbers: 21.60. Jz,21.60.Ev,24.75.+i, 25. 85.-W 



INTRODUCTION 



Interest in fission has recently increased since it is pro- 
posed to be used in new applications, such as accelera- 
tor driven systems, new electro-nuclear cycles - thorium 
based fuel cycle -, and the next generation of exotic beam 
facilities. For these applications, there is an important 
need for fission cross-sections in a large range of excita- 
tion energies, and also for mass-charge fission fragment 
distributions. For instance, precise knowledge of produc- 
tion rates of secondary long-lived fission residues and of 
neutron-rich isotopes is crucial for designing and simu- 
lating these new facilities. It is worth pointing out that 
relevant measurements of mass and charge distributions 
have been performed recently. For instance the produc- 
tion of exotic nuclei has been measured from spallation 
reactions of 1 A GeV 238 U projectiles on an hydrogen 
target pj and isotopic yields have been deduced for ele- 
ments between 58 Ni and 163 Eu. Furthermore, thanks to 
secondary beam facilities, fission properties of 70 short- 
lived radioactive nuclei can be found in references 0, 0] • 
Such a systematic analysis of the fission properties covers 
a wide region of the nuclide chart and the transition be- 
tween single- and double-humped mass distributions has 
been observed with a triple-humped structure for 227 Th. 
It is important to test the accuracy of the theoretical 
prediction using the data in order to gain confidence in 
its predictions when applied to widely extended domains 
such as fission of nuclei far from stability, and fission for 
a large range of excitation energies. 
From a theoretical point of view, the description of fis- 
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sion process stands at the crossroads of many subjects 
in the forefront of research. Both static and dynamical 
properties of the fissionning system are required, namely: 
nuclear configurations far from equilibrium, the inter- 
play of collective and intrinsic degrees of freedom, and 
the dynamics of large amplitude collective motion. The- 
oretical works generally focus on the static part of the 
fission. For instance, many studies have been devoted to 
multi-dimensional potential energy surfaces 0, from 
which fission barriers are extracted and to nuclear con- 
figurations at scission and associated fragment distribu- 
tions 0, Q • On the other hand there are very few dy- 
namical studies of fission, although dynamical effects are 
expected to play an essential role in particular in the 
descent from saddle to scission. Fragment mass distribu- 
tions have recently been obtained by solving the classical 
three-dimensional Langevin equations 8] . The influence 
of the mass asymmetry degree of freedom on the variance 
of the mass distribution has been highlighted. Two types 
of microscopic quantum dynamical calculations have also 
been performed in the past. First, in 1978 the time- 
dependent Hartree Fock method has been applied to 
fission. Second, time-dependent calculations based on 
the Generator Coordinate method using Hartree-Fock- 
Bogoliubov states have been performed, and the most 
probable fission configuration of 240 Pu has been ana- 
lyzed [Tof . The present study is an extension of this pio- 
neering work. 

In the present work we have chosen to derive the col- 
lective dynamics of fission using a time-dependent for- 
malism based on the Gaussian Overlap Approximation 
of the time-dependent Generator Coordinate (GC) the- 
ory. An alternate method could have been to first deter- 
mine the stationary solutions of the GC equations within 
the relevant domain of generator coordinates with appro- 
priate boundary conditions. The solutions of the time- 
dependent GC equations would then be expressed in a 
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straightforward manner. However, the precise form of 
the boundary conditions to be used is difficult to obtain 
when more than one generator coordinate are employed. 
Applying a time-dependent method allows one to avoid 
this problem. The only input of the calculation is the 
collective wave function chosen at t=0. Spurious reflec- 
tions of the time-dependent collective wave function on 
the edge of the finite domain are eliminated using a stan- 
dard absorption technique as explained in Section InTl 
In this paper, we focus on low energy fission-fragment dis- 
tributions of 238 U, and also on several physical aspects 
that can be clearly analyzed in this even-even fissioning 
system. Let us recall that, at low energy, elongation and 
asymmetry degrees of freedom are among the most rele- 
vant ones and that the adiabatic assumption is to a large 
extent justified |Tlj . As we explain below, time evolu- 
tion in the fission channel is described in terms of a wave 
function of Hill- Wheeler type. The latter is taken as 
a linear combination of Hartree-Fock-Bogoliubov (HFB) 
solutions characterized by the two collective degrees of 
freedom just mentioned. It is worth pointing out that 
this work relies only on the D1S effective interaction used 
at Bruyeres-le-Chatel. 

The calculation proceeds in two steps: the potential en- 
ergy surface and the collective inertia are determined 
from the first well to scission and then, the dynami- 
cal treatment of fission is performed using an approx- 
imate Time Dependent Generator Coordinate Method 
(TDGCM). Potential energy surfaces and associated 
collective inertia tensors are calculated using the con- 
strained Hartree-Fock-Bogoliubov approach with the 
D1S finite-range effective force 03 • Fission wave 
functions at time zero are constructed from the quasi- 
stationary collective states in the first well. Their time 
evolution is calculated numerically by discretizing on a 
mesh a time-dependent Schrodinger-like equation. Mass 
distributions are derived from the flux of the wave func- 
tion through scission at given A^/A^ fragmentations. 
The present work is organized as follows. The HFB for- 
malism and the TDGCM method are presented in Sec- 
tion [n] and numerical procedures are detailed in sec- 
tion lllll Section llVl is devoted to the static results, where 
the potential energy surface and pairing correlations are 
discussed. A first estimate of kinetic energy and fragment 
mass distributions, obtained from a "static" calculation 
at scission, are discussed. Mass distributions obtained 
from the full time-dependent calculations are presented 
in Section [V] and the crucial role played by dynamical 
effects is analyzed. 



acterizing the shape evolution of the nucleus. In a self- 
consistent formalism these shapes can be generated by 
means of external fields represented by the operators: 
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These moments govern mass axial deformation and left- 
right asymmetry of the nucleus, respectively. For well- 
separated fragments one can express the mean values of 
these operators in terms of < Q§) > (< Q^ Q >) and 
< Qw > (< Q30 >)) the mean quadrupole (octupolc) 
deformations of the heavy and light fragments, respec- 
tively, d m the distance between their centers of mass, 
and [i the reduced mass : 



< 



?20 > - < Qw > + < Q20 > +Wm 



< Q30 > = < Q30 > + < <?30 > 

;W "' -{A H <Q L W >-A L <Q% >) 



A H + A L 

3 A H - A L 



-2K 



Au + Ar, ' 



with 



A H A L 



Ah+Al 



(3) 



(4) 



Relations (0 and 10} have only been used in the present 
work to check the validity of the computer program for 
configurations close to scission. 

The intrinsic axially-deformed states |3>(<?20j <?3o)) of the 
fissile system are taken as the solutions of the constrained 
Hartree-Fock-Bogoliubov variational principle [L3 : 

<S($(<Z2o, q 3 o)\H-X N N-X z Z-J2 \Qi\HQ20, = , 

i 

(5) 

the Lagrange parameters Ajv, Xz, and being deduced 
from: 

($(g20,93o)|-N]$(g20,93o)) = N, 



($(920, g3o)|^|*(g20, 930)} = Z, 



(6) 



II. FORMALISM 



In low energy fission the adiabatic hypothesis seems 
to be justified [ll| and, therefore, collective and intrin- 
sic degrees of freedom can be decoupled. Furthermore, 
we assume that the collective motion of the system can 
be described in terms of a few collective variables char- 



and 



(920 , 930 ) I <9i I $ (<720 , 930 ) ) 



In Eq. JSJ), Qi is the set of external field operators (Q20, 
Q3O1 Q10), an d H is the nuclear many-body effective 
Hamiltonian built with the finite-range effective force 
D1S |l3j . The additional constraint on the dipole mass 
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operator is used in order to fix the position of the center 
of mass of the whole system. This is accomplished by 
setting < Qio > — 0, where: 



(7) 



i=l 



The system of Eqs. (|5J) and © is solved numerically for 
each set of deformations by expanding the single parti- 
cle states onto an axial harmonic oscillator (HO) basis. 
For small elongation, < 1720 < 190 b, one-center bases 
with N = 14 major shells have been considered, whereas 
for well-elongated configurations q 2 o > 190 b, two-center 
bases with N = 1 1 for each displaced HO basis have been 
used. Because calculations are performed in an even- 
even nucleus for which K = - with K the projection 
of the spin onto the symmetry axis -, the HFB nuclear 
states are even under time-reversal symmetry T . Fur- 
thermore, we restrict the Bogoliubov space by imposing 
the self-consistent symmetry 2Ti2, where n2 is the re- 
flection with respect to the xOz plane. Let us mention 
that the octupole operator breaks the parity symmetry. 
However, since PQzqP^ 1 = —Q30 and PHP = H with 
P the parity operator, constrained HFB calculations can 
be restricted to positive values of (730, and negative ones 
are obtained from ^(920, -930)) = P\4>{q 2 o, 93o))- 

The nucleus time-dependent state is defined as a linear 
combination of the basis states 10(920, 93o)} : 

!*(*)> = / / d 920 d<230 /(?20,?30,*) |0(<?2O,93O)} , (8) 



where f(q 2 o,q30,t) is a time-dependent weight function 
which is obtained by applying the variational principle: 



2 *(t)\H- ift~|#(t))dt = , (9) 



$f*(q20,qso,t) J tl 



St' 



where H is the same microscopic Hamiltonian as the 
one introduced in Eq. (JSJ. The result is the well-known 
Hill- Wheeler equation which reduces to a time-dependent 
Schrodinger equation when the GCM problem is solved 
using the Gaussian Overlap Approximation (GOA) |T(| : 

ij i ,n -f <%(920,930,*) 

Hcoii ff(<?20,930,£) = in- 7T, ■ (10) 

at 

The collective wave functions 3(920, 930, t) solutions of 
Eq. {TP} are related to the weight functions /(920, 930, t) 
through the following relation: 



7(320, 330,*) = f / dq 20 dq' 30 f(q 20 ,q^ ,t) 
^(920,930,9 2 o,9 3 o) > 



(11) 



where 1 'a (q 2 o, 930, 920, 9jso) is the square root kernel of the 
overlap kernel: 

-f(920,930,92 ,9 30 ) = (0(920, 93oM92o,9 3 o))- 



The exact form of the collective Hamiltonian H co u de- 
duced from the GOA can be found in 0, 0|. In the 
present derivation of this Hamiltonian, the widths G22, 
G23 and G33 of the gaussian overlap between differently 
deformed constrained HFB states have been assumed 
to be constant. Numerical calculation of these widths 
shows that they vary very slowly in the whole <720 _ 930 
domain considered here and that their variations can be 
neglected. With this assumption, the two-dimensional 
collective Hamiltonian reads: 

h 2 , 3 < d d 
Hcoii = — - V ~ — Bij(q 2 o,q3o)x 1-^(920,930) 



£ AVy(g 2 o,93o) , 



i,i=2 

(12) 

where ^(920,930) is the constrained HFB deformation 
energy, AV^ (920,930) are the so-called zero-point-energy 
corrections, and B t j ((720, 93o) is the inverse of the inertia 
tensor Aiij{q20j 930 ) associated with the quadrupole and 
octupole modes. In this work, we have taken for M^ , 
instead of the GCM+GOA inertia tensor, the one de- 
duced from the ATDHF theory with the Inglis-Belyaev 
approximation. The reason for this replacement is that 
the ATDHF theory appears to give a better account of 
the nuclear collective inertia than the GCM one. This 
question has been extensively discussed in the literature 
(see e.g. ^(| and references therein). 
The element (ij) of the Inglis-Belyaev inertia tensor can 
be expressed as: 



Mi 



U)-i(Jlf(-»>) 



,<■/ 



ih- 1 
hi ■ 



(13) 



fc,Z=2,3 



In Eq. i|13|) the moments of order -k are calculated as: 



M 



(-*) _ 



(0 (920 , 930 ) I Q*0 1 M^} (mH <2 jO I <K920 , 930 ) ) 



(14) 

where \[iv) are two quasi-particle states with ener- 
gies Efj, + E v built on ^(920, 93o)), and Q i0 the 
quadrupole/octupole deformation operator defined in 
Eqs. £[} and respectively. 

Let us mention that the collective Hamiltonian H co u in 
Eq. (|12|) is hermitian because: i) all inertia are real, and 
ii) -623 = i?32- 

In addition, from Eq. (JSJ, one finds that the collective 
Hamiltonian and the overlap kernel are even under the 
change of q 30 into — 930. Hence, Eq. (fTUfl propagates the 
collective wave function g{q 2 $, 930, t) without mixing par- 
ity components. In particular, if the initial wave func- 
tion g(q 2 Q, 930, t = 0) has a good parity ir, the full time- 
dependent state Eq. (fill) will be an eigenstate of P with 
eigenvalue tt. 

It is important to emphasize at this stage that the ap- 
proach presented here requires only the use of an effective 
force. We recall that the interaction D1S permits one to 
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employ the full HFB theory and consequently to treat 
the mean field and the pairing correlations on the same 
footing at each deformation. Also, the collective Hamil- 
tonian H co u , as derived from the GCM procedure, is fully 
microscopic and relies exclusively on the interaction D1S. 
Finally let us also add that the original D1S force is used, 
which means that no readjustment of the parameters has 
been made for the application reported in this paper. 
The collective Hamiltonian extracted with our proce- 
dure looks like those employed in phenomenological ap- 
proaches. However, the form used in the present work 
directly follows from the TDGCM theory and the GOA 
ansatz. We emphasize that the 2x2 inertia tensor de- 
pends on the coordinates and is non-diagonal. Since this 
situation has not been much studied, numerical meth- 
ods used to solve Eq. (|10|l are presented in Section \LU\ 
They differ from the ones previously discussed in ref. |l(j 
because of, first the large domain of deformation consid- 
ered here, and second the symmetries of the constraints, 
which lead to numerical uncertainties when implement- 
ing the previously-used procedures. 



III. NUMERICAL METHODS 

A. Discretization of the collective variables 

In order to preserve the hermiticity of the collective 
hamiltonian, the discretization of the collective variables 
has been performed by expressing the double integral of 
the functional: 

F(t) = dq 20 dq 30 g*(q 2 o,q30,t) 

J J d ( 15 ) 

(Hcoii - ifr 5(920, 930, t) , 
with finite differences: 

F{t) = J2a*{i,k,t) Kaji g(j,l,t), (16) 

ik.jl 

and by deriving the discretized equation from the varia- 
tional principle: 



dF(t) 
dg*(i,k,t) 



= 0. 



(17) 



In Eq. (TSJ, K is the symmetric matrix representing H co u 
(whose full expression is given in Appendix A), and the 
labels i, k and j, I correspond to the q 2 o & n d 930 variables 
through q 20 (i) = (i-l)*Aq 20 , and q 3 o(k) = (k-l)*Aq 3() . 



The time-dependent GCM+GOA equation becomes: 

dt 



Y,K iktjl g(j,l,t)= ih — g{i,k,t). (18) 
ji 



In practice, the two-dimensional discretized form 
Eq. (|18fl has been reduced to a one-dimensional prob- 
lem by defining a linear index m = I + (k — \).l max , with 



Imax the largest value of I on the grid, which yields: 

d 



^ K nm g m (t) = ih Q^9n{t) 



(19) 



The resulting 2x2 discretized Hamiltonian matrix H mn 
is symmetric and the hermitian character of the kinetic 
energy operator is preserved. From a numerical point of 
view, H mn is a sparse matrix. The corresponding non- 
zero elements are stored using the " row-indexed storage" 
method |l7j| . an efficient technique for reducing comput- 
ing times. 



B. Time evolution 

In matrix form, the evolution of g between t and t+At 
can be written: 



g(t + At) = e- 1 — g(t). 



(20) 



Using the Crank-Nicholson method [l7l Hfij . a unitary 
and stable algorithm, Eq. (|2*0f becomes: 



G(t + At) 



1 - i 



KAt 
2h 



1 + i 



2h 



G{t) + 0({KAtf) . (21) 



This equation can be transformed into the linear system: 

KAt. , A . KAt, , , 

(i + ^^H^ + AiHa-*— ). 9 (i) . (22) 

In this study, Eq. (|22J) is solved by successive iterations 
until convergence. The wave function g(t + At) at time 
t + At is determined from the previously known wave 
function g(t) at time t as follows: 



g {n=a) {t + At) 
+ At) 



9(t) 
(1- 



KAt. 



KAt 

''~2h 



2h 



9(t) 



(23) 



g {n) (t + At) 



Eqs. I|23|l are solved in a 520 - 930 box of finite exten- 
sion assuming 3(920, 930) t) — along the edges of the 
box. This boundary condition leads to unphysical reflec- 
tions of the time-dependent wave function on the q 2 o = 
q20max edge of the box. In order to eliminate these un- 
physical reflections the same technique as that detailed 
in ref. |l3j | has been implemented: the wave function is 
progressively absorbed in the interior of a rectangular re- 
gion q' 20max < q 20 < 9"20maa,beyond the q 20 = q 20m ax 



edge (in the present work q 2 Q T/ 



550 b, 



920r, 



= 800 



b and q" 2 om.ax — 1300 b). Inside this region, the wave 
function g (t) is multiplied at each time-step At by the 
function of Woods Saxon structure: 



F(q 20 ) 



1 



1 + e.xp(-0.015(g 2 o - H50)) 



(24) 



As mentioned in ref. |13|| . this technique is similar to 
adding an imaginary potential —ihF(q 2 o)/ At beyond the 
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boundary <?20 = Qwmax- Since At occurs in this imagi- 
nary potential, F(q 20 ) is optimized for each time-step. 
In the present study the numerical values in Eq. Ij24(l 
have been optimized to avoid reflections for a time step 
At = 1.3 * 1CT 24 s. 

The initial wave function g(t — 0) is described in terms 
of quasi-stationnary vibrational states localized in the 
first well of the potential energy surface. The states in 
question are in fact taken as the eigenstates of a modi- 
fied two-dimensional q 20 " 930 potential, where the first 
fission barrier is extrapolated to large positive values as 
mentioned in Refs. [19|. Only the states lying between 
the top of the inner barrier and 2 MeV above have been 
considered in the present work. 

Fragment mass distributions Y(Ah) are derived by a 
time-integration of the flux J(q 2 o, 930, t).nds of the wave 
function through scission at a given fragmentation: 

Y{A H )= [ dt J{q 20 , g 30 , t).nds , (25) 
Jo 

where T is the time for which the time-dependent flux 
is stabilized along the scission line. In Eq. I|25l) . n is a 
vector normal to the scission line, and J is the current 
defined from the continuity equation: 

^15(920, 930, t)\ 2 = -div,J(q 2a ,q 30 ,t). (26) 

The current J = (J 2 , J 3 ) as calculated with the collective 
Hamiltonian defined in eq. (|12J) takes the form: 

it *\ h t * n d 9 o d 9* 

^2(920, 930,£) = rr(3 B 22 ^ 9 B 22^ — 

2« oq 20 oq 20 

, * u d 9 d d 9* \ 
+9 B 23 - 9 B 23 -^ ) , 

(27) 

h dg dg* 

^3(920, 930, t) = 7T\9'B 33 — 9B33 - — 

2i dq 3Q dq 30 

, * o d 9 n d 9* \ 
+g B 32 - gB 32 - — ) . 

0920 O920 

Expression (|27|l reveals in particular that the component 
of the current in one direction involves the gradients in 
all directions. This observation will be used in sectionlvl 
where we discuss the contributions of interference terms 
between components of different parities in the initial 
state. 



IV. STATIC RESULTS 

A. Potential energy surface 

HFB calculations have been performed for 238 U with 
constraints on both the quadrupole and octupole mo- 
ments, using the mesh sizes A.q 20 = 5 - 10 b and 
Aq 3a = 2-4 b 3 / 2 . The range of investigation extends 
from spherical shapes (q 2 o = b) up to elongations of the 



exit points, which vary from q 20 „ lax = 320 b for the most 
asymmetric fission (q 30 = 44 b 3 ' 2 ) up to q 20max — 550 b 
for symmetric fragmentation (q 3 o — b 3 / 2 ). For each 
value of the quadrupole moment, the HFB calculations 
have been restricted to solutions whose excitation ener- 
gies are at most 30 MeV above the ground-state. For 
values of q 2 o near scission, this condition leads to a 
maximum value of q 3 o = 120 b 3 / 2 . HFB solutions for 
120 < q 30 < 200 b 3 / 2 have been extrapolated. 
Fig. n shows the most significant part of the HFB poten- 
tial energy surface as a function of the quadrupole and 
octupole moments. For practical reasons, the domain of 
the plot is restricted to < q 2 o < 320 b and < 930 
< 72 b 3 / 2 and energies are truncated to 25 MeV. As ex- 
pected in this actinide nucleus, the ground-state is found 
to be deformed with q 2 o « 30 b and a super-deformed 
minimum appears for an elongation close to q 2 Q = 80 b. 
Beyond this second well, two valleys appear. They are 
separated by a ridge for well-elongated shapes and lead 
either to the symmetric or to the most probable asym- 
metric fragmentations. 




FIG. 1: HFB potential energy surface as a function of q2o and 
q30 collective variables in 238 U. 



For each asymmetry, the determination of scission con- 
figurations is made by increasing the elongation step by 
step: the constrained HFB wave function at a given q 2 o 
is generated from a previous solution at a slightly lower 
elongation while keeping q 3 Q fixed. This method relies 
on the scission mechanism studied in [10] . It is assumed 
that scission occurs for a given value of g3o when the sys- 
tem falls from the so-called " fission valley" to the " fusion 
valley" describing well-separated fragments. 

The main criterion used to define exit points and to 
separate pre- and post- scission configurations is ob- 
tained by looking at the nucleon density in the neck: 
we consider that the system is composed of two frag- 
ments when the density in the neck is less than 0.01 
nucleon/fm 3 . This is illustrated in Fig. [5J where den- 
sity contours are plotted for a given asymmetry q 3a = 
44 6 3 / 2 and an increasing elongation. Contour lines are 
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separated by 0.01 nucleon/fm 3 . Figures Ufa) and H£b) 
correspond to pre-scission configurations and figure [^c) 
to a post-scission one. Let us note that the two criteria 
described in Rcf. are also satisfied: a ~ 15 MeV drop 
in the energy of the total system, and a ~ 30% decrease 
of the hexadecapole moment are observed when scission 
occurs. 
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FIG. 2: Proton plus neutron density contours at a given asym- 
metry <73o = 44 b 3//2 for different elongations a) g2o = 310 b, 
b) 920 = 320 b and c) 920 = 330 b. Contour fines are separated 
by 0.01 nucleons fm -3 . 



It is worth pointing out that the constrained HFB 
method does not impose an a priori shape to the fission- 
ing system. All types of deformations which are not im- 
posed take the values that minimize the total nuclear en- 
ergy with both the nuclear mean field and pairing field de- 
termined self-consistently. Results concerning fragment 
deformations at scission will be presented in a forthcom- 
ing publication. 

Near the exit points, the z-location of the neck, z nec k, 
is determined as the z-value for which the nucleon den- 
sity integrated over r is minimum. Properties of the 
fragments, such as their masses and their charges, their 
deformations and the distance between their centers of 
charge are calculated from integrations in the left and 
right half-spaces on either sides of the z = z ne( ± plane. 
As an example, the distance d between the centers of 
charge of the fragments is plotted in Fig. [3] as a function 
of the heavy fragment mass. It is found to be maxi- 
mum for Ah = 119 with d = 20.27 fm and minimum for 
Ah = 134 with d = 15.88 fm. Precise values of this frag- 
ment center of charge distance are crucial because they 
govern the Total Kinetic Energy (TKE) distribution, as 
discussed in Section IIV CI As a test, we have checked 
that the analytical relations in Eq. © are fulfilled. 



E 20 




120 125 130 135 140 145 150 155 160 

Heavy fragment mass 



FIG. 3: Distance between the centers of charge of the frag- 
ments as a function of the heavy fragment mass. 



B. Pairing correlations 

Fig. ^] shows the pairing energy E pair = |Tr(A«;), 
where A and k are the pairing field and the pairing ten- 
sor, respectively. We clearly see that pairing is not con- 
stant as a function of elongation. As expected, minima 
are found inside the wells and maxima at the top of bar- 
riers. Furthermore, the total pairing energy, E pa i r , is 
predicted to be lowest in the asymmetric valley (E pa ir ~ 
6 MeV) and much larger in the symmetric one {E pa ir 
> 15MeV). These variations of the pairing correlations 
are very important since they strongly influence both the 
collective flux and the occurrence of intrinsic excitations, 
as is now explained. 




50 100 150 200 250 300 350 400 450 500 550 

q 2 „ (b) 

FIG. 4: Pairing energy as a function of q2o along the asym- 
metric (solid line) and the symmetric (dashed line) fission 
paths in 238 U. 



First, the collective inertia is known to be very sensi- 
tive to pairing correlations. The three components B22, 
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-B33 an d B23 of the inertia tensor in Eq. (|12f> are plot- 
ted in Figs.|£ta),EKb) andOJc), respectively, as functions 
of the elongation along the symmetric (dotted line) and 
asymmetric (solid line) paths. The two components B22 
and .B33 are found to be larger in the symmetric valley 
than in the asymmetric one (up to a factor of two at large 
elongation). Furthermore, whereas the non-diagonal in- 
ertia component B23 is zero for q^o = by definition, 
B23 is found to be non negligible as soon as the system 
spreads widely in the asymmetric valley. The coupling 
brought by B23 between the ^20 an d 930 modes indi- 
cates that, as time evolves, the two collective degrees of 
freedom exchange energy, which will affect, among other 
things, the kinetics of the fission process. 




600- a) 

400 - 
200 - 



50 100 150 200 250 300 350 400 450 500 550 

q 3 „ (b 3 ' 2 ) 

FIG. 5: Components of the inertia tensor a) B22 (in 
MeV b 2 h 2 ) b) B 3i (in MeV b 3 h 2 ) c) B 2S (in MeV b 5/2 ft 2 ) 
as functions of 520 along the asymmetric (solid line) and the 
symmetric (dashed line) fission paths 



Second, pairing correlations characterize the amount 
of superfluidity of the collective flux and the onset of dis- 
sipation, in particular between the saddle point and the 
exit point. In the HFB approach, dissipation requires the 
creation of two quasiparticle excitations, that is a trans- 
fer of energy from the collective motion at least equal to 
2A, where A is the energy necessary to break a correlated 
pair. One expects that small values of A will favor "dis- 
sipation" . However, the excitation of the intrinsic struc- 
ture also depends on the coupling between collective and 
intrinsic degrees of freedom, which is largely unknown. 
For this reason, the question of dissipation effects will 
be addressed in future work. The proton and neutron 
gaps 2A p , 2A„ are plotted in Fig.[fJ]as functions of elon- 
gation along the asymmetric path. The corresponding 
potential energy curve is also plotted (dotted curve) to 
guide the eyes. For proton pairing correlations, we find 
2A p = 2.3 MeV at the top of the second barrier. This 
value app ears to be in good agreement with experimental 
data |20l l2ll |22]| . As a matter of fact, manifestations of 
proton pair breaking are observed in 238 U and 239 U nu- 
clei for an excitation energy of 2.3 MeV above the barrier: 



first the proton odd-even effect observed in the fragment 
mass distributions decreases exponentially for excitation 
energy slightly higher than 2.3 MeV [2(| and second, the 
total kinetic energy drops suddenly |2ll l22| . 




40 80 120 160 200 240 280 320 

q 2 „ (b) 

FIG. 6: Twice the proton (solid curve) and neutron (dashed 
curve) lowest quasi particle energies and potential energy 
(dotted curve) along the asymmetric path as functions of the 
elongation in barns. 

In Fig. we also see that the proton gap decreases 
rapidly during the first part of the descent beyond the 
saddle point, for instance 2A p = 1.4 MeV for 520 = 180 b. 
However, experimental facts show that increasing the ex- 
citation energy from to 2.2 MeV above the barrier does 
not modify the proton odd-even effect. From our point 
of view, this could be an indication that no proton pairs 
are broken during the descent from saddle to scission in 
low energy fission for excitation energies below 2.2 MeV 
above the barrier. Our interpretation is that the excita- 
tion energy supplied during the descent is shared among 
the collective degrees of freedom and not among intrin- 
sic excitations. This experimental observation gives us 
some confidence that the neglect of the coupling between 
collective and intrinsic degrees of freedom is a reasonable 
approximation to start with in low energy fission. 

Finally, no strong odd-even neutron effects are ob- 
served for the fragment mass distributions measured 
in the photofission of 238 U, regardless of excitation en- 
ergy |21|. In our calculations the neutron pairing gap 
is found to be much lower than the proton one, except 
for 160 < q 2 o < 190 b , as displayed in Fig. HO At the 
top of the second barrier the neutron gap is only 2A„ 
= 1.6 MeV. This tends to indicate that neutron pairs 
are more likely to be broken than proton ones in the 
even-even 238 U nucleus at low excitation energy. But no 
definite comparison with experimental data can be made 
since a precise knowledge of the neutron number of the 
fission fragments is made extremely difficult by the neu- 
tron evaporation. All these remarks concerning pairing 
correlations a posteriori illustrate the fact that pairing 
correlations play an essential role and that they should 
be introduced in dynamical studies of fission. 
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C. Total Kinetic Energy distribution 



D. "One-dimensional" fragment mass distribution 



As a first estimate, the total kinetic energy TKE of 
the fragments can be roughly calculated as the Coulomb 
potential energy TKE = s d Le , with d the distance 
between the centers of charge of the fragments at scis- 
sion. 

Theoretical values calculated along the scission line are 
shown in Fig. [7| as a function of Ah, the heavy frag- 
ment mass. They are compared to experimental data 
obtained from the photofission of 238 U using 6.2 MeV 
bremsstrahlung 7 rays, corresponding to an excitation 
energy close to the inner fission barrier height |2l|. We 
first notice that the general trend of the distribution is 
rather well reproduced, with a dip at Ah = 119 and a 
peak for Ah = 134. Symmetric and asymmetric wings 
are surprisingly close to experimental data. The agree- 
ment indicates that our microscopic approach, together 
with the prescription explained in Section flV Al is able to 
give a realistic description of scission configurations. The 
main difference with experimental data occurs in the re- 
gion of the most probable asymmetric fission where the 
theoretical results overestimate TKE values by ~ 6%. 
This discrepancy mainly comes from the fact that the 
nuclear contribution entering the mutual energy between 
the two fragments is not strictly zero for the correspond- 
ing scission configurations. Furthermore, the - attractive 
- exchange Coulomb energy between the fragments has 
been neglected. These two effects could lead to a decrease 
of TKE values that may reach 10-15 MeV. 




120 125 130 135 140 145 150 155 160 
Heavy fragment mass 



FIG. 7: Total kinetic energy distributions as functions of the 
heavy fragment mass. Dots indicate experimental data |21fj 
and the continuous line present predictions. 



As a first approximation, mass distributions can be de- 
rived using the fragmentation model detailed in ref. |28j | . 
Namely, collective stationary vibrations along the sole 
mass-asymmetry degree of freedom for nuclear configu- 
rations just before scission are studied. The probability 
of occurrence of a mass asymmetry (Ah, Al) correspond- 
ing to a value q 3 $ of the octupole moment is then taken 
as: 



Y(A H ,A L ) = \*+\q 30 ) 



(28) 



where "J J 1 is the positive parity eigenstate with lowest 
energy of the one-dimension collective Hamiltonian H' co u 
in the q 3 o variable: 



' coll (<?30 ; 



a 

dq. 



30 



) 
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1 



d 



2 dq 30 M 3 (q 3 o) dq 30 
-V(qao) - AV 3 (q 30 ) . 



(29) 



Here, V(q 3 o) is the HFB deformation energy along the 
scission line q 2 o = 5lo = filso), M 3 (q 3 o) is the collec- 
tive inertia, and AV3(g3o) the zero-point-energy correc- 
tion (ZPE) . The Hamiltonian of Eq. is derived from 
the usual GOA reduction of the one-dimensional Hill- 
Wheeler stationary equation obtained by taking for the 
generator coordinate the curvilinear abscissa s(q 2 o,q 3 o) 
along the scission line. A change of variable is then per- 
formed in order to express all quantities as functions of 
q 3 o. It is easy to show that the inertia and ZPE correc- 
tion appearing in Eq. (|29|l are related to those entering 
the full two-dimensional Hamiltonian 11211 : 



df 



df 



M 3 (q 30 ) = {-^—Y M 22 + 2(-i-) M 23 + Al 



dq, 



30 



dq, 



33 - 



30 



A^sfeo) = 



£3(430) 

2M 3 {q 3 vY 



with 



£3(930) = ( 



dq, 



■G 



22 



30 



2(^)G 23 
a<730 



G 



33 - 



(30) 



(31) 



(32) 



where My are the inertia defined in Eq. (|13|) and Gy 
the components of the overlap tensor calculated in the 
cranking approximation using the moments of Eq. (I14|l . 

Clearly, the model based on the Hamiltonian 129|) 
amounts to ignore all the effects of the dynamics along 
the elongation degree of freedom from the first well to 
scission. We call the mass distribution obtained in this 
way a "one-dimensional" mass distribution. 

The HFB potential energy V(q 3 $) calculated along the 
scission line is plotted in Fig. [S] as a function of the oc- 
tupole moment. The lowest energy is obtained for 530 = 
± 44 b, corresponding to the most probable fission. A 



9 




-160 -120 



FIG. 8: Potential energy along the scission line as a function 
of the octupole moment. 



secondary minimum is found for q^Q = b. These two 
wells are separated by a 11 MeV high barrier. 

"One-dimensional" distributions are shown in Fig. 
where the mass yield Eq. (|28|) is plotted (solid line) to- 
gether with the Wahl evaluation (dashed line) for 46 keV 
neutron induced fission on 237 U [24J . The maxima of the 
theoretical curve occur at Ah = 134, Al = 94 values cor- 
responding to the minima of the potential energy along 
the scission line. The fact that the experimental curve 
maxima lie close to these values indicate that the most 
probable fragmentation is due essentially to the proper- 
ties of the potential energy surface at scission, i.e mainly 
to shell effects in the nascent fragments. However, the 
"one-dimensional" approach does not reproduce neither 
the experimental peak-to- valley ratio nor the experimen- 
tal widths of the distributions - the theoretical widths 
are twice smaller than the Wahl evaluated ones-. 



Eq. H29fl with lowest energy 'F^ 1 has so far been con- 
sidered. This is certainly an oversimplifying assumption 
since, the wave function describing the collective evolu- 
tion of the nucleus will undoubtedly possess a more com- 
plicated structure at the time of reaching the scission 
line. In particular, as mentioned in ref. [2fij |. n > states 
solutions of Eq. 

may become excited due to 
the interaction between (720 and g3n degrees of freedom. 
Let us mention that, in Ref (2^ the population of the 
eigenstates has been assumed to follow a Boltzmann law 
governed by a temperature parameter. In Ref. j2(J, the 
elongation degree of freedom has been introduced using 
a classical approximation. 

Amplitudes of the first six collective states <F£, 
n = ... 5, solutions of Eq. are displayed in Fig. ITUI 
as functions of the fragment mass. As is well-known, 
such *F^ states are eigenstates of the parity operator II 
with eigenvalues n. Positive and negative parity states 
are plotted in solid and dotted lines, respectively. Each 
pair of 7r = +1 and tt — -1 levels is degenerate in en- 
ergy because the potential is symmetric with respect to 
the 530 — ► -930 transformation, and because the barrier 
between the two asymmetric wells is high (11 MeV) (see 
Fig. |SJ • We observe in Fig. ^3 that the excited states 
are more spread over mass than the ground state. For 
example, the wave functions 'Fj 1 and 'F^ -1 displayed 
in Fig. 110b ) display non-zero values up to Ah ~ 156 
whereas the lowest energy wave functions "F^ 1 and "F^ 1 
in Fig. llOfa ,) are localized in the domain 132 < Ah < 144. 
Therefore, we can expect that the introduction of these 
excited states in the definition of the mass yield Eq. I|28() 
will broaden the mass distribution in the asymmetric re- 
gion. 




60 72 84 96 108 120 132 144 156 168 180 

Fragment Mass 

FIG. 9: Comparison between the "one-dimensional" mass 
fragment distribution obtained from 1281 (solid line), and the 
Wahl evaluation (dashed line) |24| . 
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FIG. 10: Amplitudes of the first six collective states in the 
asymmetry variable at scission as functions of the fragment 
mass. Positive (negative) parity states are plotted in solid 
(dotted) lines. Excitation energies, measured from the low- 
est energy state, are a) E = MeV, b) E = 3.09 MeV and, 
c) E = 4.76 MeV. 



One must note however that only the solution of 
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V. DYNAMICAL RESULTS 
A. Initial states 

The time-dependent evolution of the system has been 
calculated using different initial conditions in the first 
well of the potential energy surface. Calculations have 
been performed from t = up to maximum times for 
which the flux of the time-dependent collective wave 
function along the scission line has become stabilized. 

We first discuss the effect of the structure of the ini- 
tial state on the mass distribution. In order to define 
the initial conditions we imagine that the nucleus is 
a compound system described in terms of complicated 
quasi-stationnary states which decay into various chan- 
nels (neutron and 7-ray emission and fission). In the 
case of the even-even K = 238 U nucleus studied here, 
we assume that states which decay through fission can 
be described by the simple form: 

|*p,K=o,/,m) = (27r)-4F / , M (0) J dq fZ(q,t= 0) |0(q)}, 

(33) 

where Yj M (fl) are spherical harmonics and SI the Euler 
angles relating the intrinsic axes of the nucleus to the 
laboratory frame of reference. 

The parity quantum number P is related to the intrin- 
sic parity tt by the following relation : 

P = n(-lY, (34) 

where I is the spin of the fissioning system. In Eq. (13311 , 
q is the set of all relevant nuclear collective deformations 
which, in the present work (see Eq. ©) is restricted to 
(920,930)- 

As already mentioned in Section ^ initial states 
fl£(<?20, 930 ,t = 0) (related to the f%(q 20 , 930, t = 0) func- 
tions as in Eq. JTT}) are taken as eigenstates of the modi- 
fied two-dimensional first well V'(q 2 Q, 930)) where the po- 
tential has been extrapolated at large deformations as 
shown in Fig. ^2 They are solutions of the equation: 

H'aoii ^(920, 930, t = 0) = E* ^(920, 930, t = 0) , (35) 

where H' co u is the Hamiltonian defined in Eq. 112(1 with 
^(920, 930) replaced by V'(q 20 , 930)- 

Because H' cM (q 20 , 9so , gf^ , gf^ ) 

HcoiA&o, -930, gf^, -gf^) tnese initial states can 
be chosen as eigenstates of the parity operator with 
eigenvalues tt = ± 1 p7| : 

fl£(920,-930,* = 0) = tt , 9^(920,930, t = 0) . (36) 

The potential curves V(q 2 o, 930 = 0) and V'(q 2 o, q^a = 0) 
including zero-point-energy corrections are displayed in 
Fig. El The eigen-energies of H' cM are also shown. Ex- 
citation energies of the compound nucleus in the interval 
[ Bj , Bj +2.5 MeV ] will be considered in this work 



focussing on low-energy fission, where B\ is the first bar- 
rier height. In this energy range, the mean-level spacing 
is found to be around 130 keV. Therefore, 19 states are 
possible initial candidates for our dynamical calculations. 
All these states are located above the outer symmetric 
saddle point but below the outer asymmetric one. They 
correspond to multi- quadrupole and octupole phonons, 
and have different components along the (730 and 520 di- 
rections. Significant effects on the fragment mass distri- 
butions are mainly due to the parity of the initial states. 




10 20 30 40 50 60 70 80 90 100 110 120 



q 2 » (b) 

FIG. 11: Potential curves F(q20,930 = 0) and V'(g2o, 930 = 0) 
including zero-point-energy corrections (continuous and dot- 
ted lines, respectively) , and collective eigenstates of the mod- 
ified F'(g20,930 = 0) potential (horizontal segments). 

In Fig. ^1 fragment mass distributions, calculated 
with formula (|25H and initial states of definite parity, are 
plotted separately. The solid and dotted curves corre- 
spond to initial states whose intrinsic parity is positive 
or negative exclusively. They are located at 2.4 MeV and 
2.3 MeV above the first barrier, respectively. 



10V 




Fragment mass 

FIG. 12: Fragment mass distributions obtained for initial 
states having positive parity (solid line) and negative parity 
(dotted line). 

We see that the main difference between the two re- 
sults is the peak-to-valley ratio which is of the order of 



11 



50 for the positive parity initial state and infinite for the 
negative one. The fact that no symmetric fission is found 
when the initial state has a negative parity is due to the 
fact that <7^ 1 (g20 J 930 = 0, t) = if tt = -1 for any time t. 
As a consequence, the flux of the wave function through 
the scission line at (730 = vanishes. 
In the applications presented below we use initial states 
that don't have a definite parity. As we observed at the 
end of section ITTT1 there are interferences between states 
of different parities in the calculation of the flux. How- 
ever, due to the symmetries of the inertia tensor, these 
interferences do not contribute to the symmetric mass 
fragmentation. ^From the previous discussion we infer 
that our predictions of symmetric fission will be affected 
by the proportions of collective states with negative and 
postive intrinsic parity. In order to get an estimate of 
these proportions we assume that they are the same as 
in the compound sytem n+ 237 U. More precisely, by us- 
ing Eq. I|34l) we define fission cross sections, a(n = +1,E) 
and <j(it = -1,E) corresponding to components of intrinsic 
parity in the compound system through the relations: 

a(n = -l,E)= a C N(P,I,E)Pf{P,I,E) 

I=2p,P=-l 

+ Yl <7CN(P,I,E)P f (P,I,E) , 
I=2p+l,P=+l 

a(n = +l,E) = <tcn(PJ,E) P f (P,I,E) 

I=2p,P= + l 

+ vcN(P,I,E)P f (P,I,E) , 

I=2p+l,P=— 1 

(37) 

where E, P, <Jcn(P, I, E) and Pf(P, I, E) are the energy, 
the parity (defined in Eq. H34[l ). the formation cross sec- 
tion and the fission probability of the compound nucleus, 
respectively. 

Formation cross sections have been calculated using 
the Hauser-Feschbach theory with the optical potential 
model of Ref. [2^ and fission probabilities have been de- 
duced from a statistical model calculation p9l |. 

Then, we define probabilities by the following fractions: 



p-(E) = 



p+(E) 



ct(tt = -1,E) 



a(% = -l,E)+a(% = +l,E) ' 

g(7T = +1,E) 

a(Tr = -l,E)+a(7T = +l,E) ' 



(38) 



They represent the population of states in the com- 
pound system which have a given intrinsic parity and 
which decay to fission. It is with the help of these prob- 
abilities that we determine the mixing of parities in the 
initial states. Numerical values for the reaction n+ 237 U 
are given in Table [0 for two excitation energies as mea- 
sured from the top of the first barrier. 

Mass distributions obtained for these two energies are 
displayed in Fig. ^| the solid and dashed lines corre- 
spond to the two energies, 2.4 MeV and 1.1 MeV respec- 
tively. 



E (MeV) 


1.1 


2.4 


P+(E) % 


77 


54 


P~(E) % 


23 


46 



TABLE I: Percentages of positive and negative intrinsic parity 
states populated in the compound nucleus 238 U by the n+ 237 U 
reaction for two excitation energies. 
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Fragment Mass 

FIG. 13: Fragment mass distributions obtained for two initial 
states. Solid line: E = 2.4 MeV with p + =54 % and p~ =46 %. 
Dashed line: E = 1.1 MeV with p+ = 77 % and p~ = 23 %. 



One observes that the symmetric fragmentation is 
slightly higher at excitation energy 1.1 MeV than at 2.4 
MeV. Clearly, this approach does not reproduce an es- 
sential feature of measured or evaluated mass fragment 
distributions namely, a sensitive increase of the symmet- 
ric fission yield with increasing neutron energy. As can be 
inferred from previous discussions, this discrepancy is a 
direct consequence of the rapid decrease with increasing 
energy of positive parity components in the initial state. 
More detailed comparisons of the theoretical predictions 
with Wahl evaluations [24| Fig.ll4lindicate. however, that 
the microscopic approach reproduces successfully various 
characteristics of the mass distribution. 

For instance, the comparison at 2.4 MeV shows that 
the main features of Wahl's distribution, position and 
height of the maxima, ratio peak to valley and broaden- 
ing of the distribution as well, are satisfactorily repro- 
duced by the theory. The agreement at lower energy is 
not as good, due essentially to the discrepancy mentioned 
above. 

Initial conditions appear to be crucial for the predic- 
tion of mass distributions at low energy. In view of the 
quality of the results presented above we consider study- 
ing more carefully this question in future works. 
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FIG. 14: Theoretical mass distributions (solid lines) are com- 
pared with the Wahl evaluations of neutron induced fission of 
238 U [2^| (dashed lines). Excitation energies of the compound 
238 U nucleus measured above the barrier are a) E = 2.4 MeV, 
b) E = 1.1 MeV 



B. Dynamical effects 

In order to analyze the influence of dynamical effects, 
the fragment mass distribution obtained for the initial 
state located 2.4 MeV above the first barrier shown in 
Fig. 114b ) is compared in Fig. 1 151 with our previous "one- 
dimensional" distribution. Fig. ED also shows the evalu- 
ated data from the Wahl systematics (dashed curve) [24[ • 
We first note that the maxima of the two theoretical 
distributions are both located around Ah = 134 and Al 
= 104, in good agreement with evaluated data. As al- 

this is a confirmation 



ready mentioned in Section lIV Dl 
that the most probable fragmentation is due essentially to 
shell effects in the nascent fragments and not to dynam- 
ical effects. The widths of the peaks obtained from the 
full dynamical calculation are much larger- about twice 
as large - than those of the "one-dimensional" one, and 
consequently they are in much better agreement with the 
Wahl evaluated data. Clearly the dynamics is found to 
play a major role in the broadening of the fragment mass 
distributions. In the present dynamical calculation the 



FIG. 15: Comparison between the "one-dimensional" mass 
distribution of Fig. [5] (dotted line), the mass distribution re- 
sulting from the dynamical calculation (solid line) with the 
initial state state located 2.4 MeV above the barrier, and the 
Wahl evaluation (dashed line) |24|. 



broadening is clearly due to the interaction between the 
elongation and the asymmetry degrees of freedom, which 
results from both the potential energy and the inertia 
variations. As already discussed in Section ITTT1 these ef- 
fects are especially important in the descent from saddle 
to scission, where the inertia component B23 is found 
to be large (see Fig. EJ. In fact the cross term in the 
kinetic energy of the collective Hamiltonian appears to 
be responsible for exchanges of energy between the two 
modes and for the spreading of the time-dependent wave 
function in the asymmetric valley. 

In order to quantitatively analyze those effects, the time 
dependent wave function g(q2o, <l30,t) has been expanded 
over the "one-dimensional" states ^(^n) described in 
Section CxEl along the line q 2 o = (920)5 = /(<Z3o): 

ff((«£»).,«30,*) = K(l3o)- (39) 

The weight coefficients C n (t) can be calculated as: 

C n (t) = J dq 30 g((q 20 ) s ,q 3 o,t)K(Q3o), (40) 

and the fraction of each "one-dimensional" state con- 
tained in the dynamical solution at scission is given by: 



Rn(t) 



\c n (t)f 



£JCm(t)| 2 ' 



(41) 



Results for R n (t = 0.96 10" 20 s) are listed in table H2 in 
the case of the dynamical wave function corresponding 
to the initial state considered here. 

These results indicate that the dynamical wave func- 
tion is spread over many "one-dimensional" states ^ 



n 


1-2 


3-4 


5-6 


7-8 


9-10 


R n % 


35.2 


8.6 


36.7 


6.9 


12.6 



TABLE II: Percentages of the "one-dimensional" states 
contained in the dynamical solution at scission for 
t = 0.96 10~ 20 s. 

and that the relative contribution of the two low- 
energy states is only 35.2 %. As it appears, the "one- 
dimensional" definition l|28|) of the fragment yield is not 
pertinent, and dynamical effects should be fully taken 
into account in order to obtain realistic predictions for 
fragment mass distributions. 

VI. CONCLUSION 

In this work, we have presented a theoretical frame- 
work and numerical techniques allowing one to describe 
fission mass distribution in a completely microscopic way. 
The method is based on a HFB description of the in- 
ternal structure of the fissioning system. The collective 
dynamics is derived from a time-dependent quantum- 
mechanical formalism where the wave- function of the sys- 
tem is of GCM form. A reduction of the GCM equation 
to a Schrodinger equation is made by means of usual tech- 
niques based on the Gaussian Overlap Approximation. 
Such an approach has the advantage of describing the 
evolution of heavy nuclei toward fission in a completely 
quantum-mechanical fashion and without phenomenolog- 
ical parameters. 

Properties of the fissioning system which have a large 
influence on collective dynamics have been discussed. 
Among them, the most important is the variation with 
deformation of the nuclear superfluidity induced by pair- 
ing correlations. In addition to strongly influence the 
magnitude of the collective inertia, these correlations are 
essential in our approach because they validate the adia- 
batic hypothesis as a first approximation for the descrip- 
tion of low energy fission. 

In the present application of this method to 238 U fis- 
sion, two kinds of observables have been examined and 
compared to experimental data: the kinetic energy dis- 
tribution and mass distribution of fission fragments. The 
kinetic energy distribution, which has been derived from 
the mutual Coulomb energy of the fragments at scission, 
is found to be in good agreement with data. A small dis- 
crepancy (6%) is found around the most probable frag- 
mentation region, which could originate from the fact 
that the nuclear contribution entering the mutual energy 
between the two fragments is not strictly zero for the 
corresponding scission configurations and that the attrac- 
tive exchange Coulomb energy between the fragments has 
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been neglected. Concerning fragment mass distributions, 
the main result of the present study is that dynamical ef- 
fects taking place all along the evolution of the nucleus 
are essential in order to obtain widths in agreement with 
experimental data. In contrast, the maxima of the dis- 
tributions are determined by the static properties of the 
potential energy surface in the scission region, that is by 
shell effects in the nascent fragments. Finally, the influ- 
ence of the choice of the initial state has been studied. 
In particular, symmetric fragment yields are found to 
be strongly influenced by the parity composition of the 
initial state. The quality of the results reported here en- 
courages us to pursue further studies of fission along these 
lines, with some additional improvements. For instance, 
as suggested by the work cited in Ref . [3(j , we cannot ex- 
clude that several valleys due to other collective modes, 
such as hexadecapole or higher multipole deformation, 
can appear in some fissioning systems. Extensions to mi- 
croscopic calculations involving three or more collective 
coordinates are envisaged. 
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APPENDIX A: HAMILTONIAN MATRIX 

Starting from the functional Eq. 1151) . the matrix el- 
ements of the Hamiltonian matrix K can be expressed 
from: 

(Hcoiig) (», k, t) = J2 K *Jt 9(3, 1, t). (Al) 
jl 

By writing: 

3 

H coH = f iJ + ^(920, 930) , (A2) 

with 

h 2 d d 

3 

^feo,93o) = ^(920,930)- ^2 AV l j(q 2 Q,q3o) , 

(A3) 
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the different terms contributing to Eq. lAlfl are: 



(A4) 



(A5) 



(T22 + f 33 )g(i, k, t) = -^-{[-B 22 (i -l,k)- B 22 (i, k)]g{i - 1, k, t) 

+[B 22 (i — l,k) + 2B 22 (i, k) + B 22 (i + 1, k)]g(i, k, t) 
+ [-B 22 (i, k) - B 22 (i + 1, k)]g(i + 1, k, t)} 
+T^^~ B ^(h k - 1) - B 33 (i, k)]g(i, k — l,t) 

+ [^33(1, k - 1) + 2B 33 {i, k) + B 33 {i, k + l)]g(i, k, t) 
+[-B 33 (i, k) - B 33 (i, k + l)]g(i, k+1, t)} , 

(f 23 + f 32 )g(i,k,t) = 1 {[B 23 (i-l,k+l)+B 23 (i,k+l) 

16Aq 2a Aq 3a 

+B 23 {i -l,k) + B 23 (i, k)]g(i - 1, k + l,t) 
+ [B 23 (i -l,k-l) + B 23 (i + 1, k + 1) 
-B 23 {i - 1, k + 1) - B 23 {i + l,k- l)]g(i, k, t) 
+ [-B 23 (i-l,k)-B 23 (i,k) 
-B 23 (i - 1, k - 1) - S 2 3( l; fc - l)]g(i - 1, fc - 1, t) 
+ [-B 23 {i, k+1) — B 23 (i + 1, k + 1) 
--B 23 (i, fc) - B 23 (i + 1, k)]g(i + l,k+l,t) 
+ [B 23 (i,k) + B 23 (i+l,k) 

+B 23 (i, k - 1) + B 23 {i + 1, k - l)]g(i + 1, k - 1, t)}, 

and 

V-(i, fe) g(i, k, t) = {V(i, k) - AV 22 (i, k) - AV 22 (i, k) - 2AV 23 (i, k)}g(i, k, t). (A6) 

The labels i and k are related to q 2 Q and 930, respectively, and Aq 2 $ and Aq3o are the associated discretization steps. 
The following approximation has been used for the inertia term: 

B r] (i ± 1 fc) » 1 (i , fc) + j (i ± 1 , k) ) , ( A7) 

and, for products -F of two functions F\ and F 2 the following prescription has been assumed: 

F(i+l,k) =F 1 (i+l,k)F 2 (i+l,k) 

2 I 2 2 (A8) 

w + 1, fc) + fc))(F 2 (t + 1, fc) + F 2 (z, fc)) . 

I 
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